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Abstract 

Data acquired by the Galileo magnetometer on five passes by Ganymede have been used 
to characterize Ganymede’s internal magnetic moments. Three of the five passes were 
useful for determination of the internal moments through quadrupole order. Models 
representing the internal field as the sum of dipole and quadrupole terms or as the sum of 
a permanent dipole field upon which is superimposed an induced magnetic dipole driven 
by the time varying component of the externally imposed magnetic field of Jupiter’s 
magnetosphere give equally satisfactory fits to the data. The permanent dipole moment 
has an equatorial field magnitude 719 nT and is tilted by 176° from the spin axis with the 
pole in the southern hemisphere rotated by 24° from the Jupiter-facing meridian plane 
towards the trailing hemisphere. The data are consistent with an inductive response of a 
good electrical conductor of radius approximately 1 Ganymede radius. Although the data 
do not enable us to establish the presence of an inductive response beyond doubt, we 
favor the inductive response model because it gives a good fit to the data using only 4 
parameters to describe the internal sources of fields, whereas the equally good dipole plus 
quadrupole fit requires 8 parameters. An inductive response is consistent with a buried 
conducting shell, probably liquid water with dissolved electrolytes, somewhere in the 
first few hundred km below Ganymede’s surface. The depth at which the ocean is buried 
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beneath the surface is somewhat uncertain, but our favored model suggests a depth of 
order 150 km. As both temperature and pressure increase with depth and the melting 
temperature of pure ice decreases to a minimum at -170 km depth, it seems possible that 
near this location, a layer of water would be sandwiched between layers of ice. 
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Introduction 

Ganymede, Jupiter’s third Galilean satellite, has an internal magnetic dipole (Kivelson 
et al. 1996, 1997) strong enough to create its own mini-magnetosphere inside of Jupiter’s 
larger one (Gumett et al. 1996, Kivelson et al. 1996). Magnetopause crossings were 
identified in the Galileo magnetometer data (Kivelson et al. 1992) acquired during flybys 
of Ganymede (Kivelson et al. 1997). Williams et al. (1997a) used data from the 
Energetic Particle Detector (EPD) (Williams et al. 1992) to confirm that energetic 
electrons are trapped on closed field lines (connected to Ganymede at both ends) within 
Ganymede’s magnetosphere. This important feature of the near-Ganymede environment 
is consistent with analysis of field line resonances found in the magnetometer data by 


Volwerk et al. (1999). 

The magnetometer data acquired near closest approach to Ganymede on the G 1 (June 
27, 1996) and G2 (September 6, 1996) flybys were well modeled as the field of a 
centered internal dipole with equatorial surface field strength approximately 750 nT tilted 
by -170° with respect to Ganymede’s rotational axis (Kivelson et al. 1997). 

With data now available from several additional passes (see Table I which lists all 
Ganymede passes including one in late 2000 and Figure 1 which shows the trajectories of 
the different passes), it is appropriate to revisit the modeling of the internal field. We 
present here several new results. First, based on an improved analysis of the data from 
selected passes, we place tighter constraints on the dipole moment and discuss limits on 
the quadrupole contributions. The quadrupole moment can be taken as evidence for the 
depth at which the dynamo driving the field is located (Lowes 1974, Elphic and Russell 
1978, Connemey, 1993), although the argument is not without its critics. We also 
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consider the possibilities of remanent ferromagnetism (Crary and Bagenal, 1998) or 
magnetoconvection (Sarson et al. 1997) as the field generation mechanism. Second, we 
address the question of whether the dipole moment changes in time or remains 
unchanged from one pass to the next. We consider the latter possibility because 
Ganymede could, in principle, respond inductively to time variations of the external 
magnetic field present at its location in the Jovian magnetosphere. Temporal variations 
arise because Jupiter’s tilted dipole moment changes its orientation as the planet rotates. 
Ganymede’s internal structure appears to include a metallic core, a rocky mantle and an 
icy outer layer, a model inferred from measurements of the gravitational moments 
(Anderson et al., 1996) and magnetic data (Schubert et al. 1996, McKinnon 1997). An 
inductive response could be present if the icy layer contains electrically conducting paths 
as, for example, in regions of partial or complete melt of sufficient thickness. 

The icy moons Europa and Callisto respond inductively to the variations of magnetic 
field present at their orbits but neither possesses a substantial permanent magnetic 
moment (Khurana et al. 1998; Kivelson et al. 1999). The measured magnetic 
perturbations observed on different passes by these bodies have been interpreted as 
evidence for subsurface oceans or analogous conducting layers (Kivelson et al., 2000). At 
Ganymede, however, extracting evidence of an inductive response in the presence of a 
large permanent magnetic moment requires a particularly refined assessment of the 
changes in the magnetic moment from pass to pass. For a perfectly conducting sphere, 
induced surface currents produce magnetic perturbations that cancel the normal 
component of a varying external field on the surface. The time-varying field at 
Ganymede is directed nearly along the Jupiter-Ganymede vector and has an amplitude of 
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~100 nT (Table 1). An induced internal dipole with 100 nT polar field (50 nT equatorial) 
antiparallel to the time-varying field along the Jupiter-Ganymede vector satisfies the 
required boundary condition at the surface. Thus the induced dipole moment will be at 
most ~50/700 or 6% of the permanent magnetic moment. As an induced magnetic 
moment is approximately perpendicular to the spin axis, induction can change the 
magnitude of the magnetic moment by no more than 0.2% and change its tilt by only 
±3.6°. 

Characterizing the internal fields of the moons of Jupiter with a high degree of 
accuracy is complicated by the fact that only a small portion of data from a flyby is 
acquired at altitudes low enough for the signature of internal sources to dominate other 
sources of magnetic field (see Table I). In addition, field perturbations arising from 
strong currents that develop within the plasma of Jupiter’s magnetosphere in the region of 
interaction with the moons are important near closest approach and must be separately 
established. 

In this paper, we first discuss our approach to fitting the data from the multiple passes 
of the Galileo mission. We then describe how we looked for an inductive response. We 
end by summarizing the implications of our findings for Ganymede’s internal structure. 

Ganymede’s internal sources fitted with dipole moments evaluated pass by pass 

The original estimates of Ganymede’s dipole moment (Kivelson et al. 1996, 1997) 
were based on data from the first two passes by Ganymede. With data now available 
from six passes by Ganymede, with one additional pass providing crucial data, we are 
able to refine the multipole moment analysis. 
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Three useful Ganymede-centered coordinate systems are illustrated in Figure 2. For 
characterizing the internal field it is conventional to use right-handed spherical 
coordinates (referred to as Gsph) that rotate with Ganymede. Longitude is measured 
from the Jupiter-facing meridian; colatitude is measured from the rotation axis. A related 
Cartesian coordinate system (£ ,t],Q has £ toward Jupiter and C, along the rotation axis. 
For characterizing the external plasma and field, it is convenient to use coordinates that 
relate to the direction of flow of the Jovian plasma. In this Cartesian coordinate system 
(referred to as GphiO), X is along the flow direction, Y is along the Ganymede-Jupiter 
vector, and Z is along the spin axis. These coordinates are analogous to the earth-centered 
GSE coordinates that relate to the direction of flow of the solar wind onto Earth’s 


Fig. 2 


environment. These definitions imply that X= -tj, Y = £ and Z = £. 

The G1 and G2 flybys occurred when Ganymede was located well above the 
magnetospheric current sheet (see Table I) in a region where Jupiter’s magnetospheric 
field points away from the planet. In order to determine whether the internal moment 
orientation changes in phase with the externally imposed field, data from passes at 
different locations relative to the current sheet must be analyzed. The G8 flyby occurred 
when Ganymede was near the center of the current sheet and the G7 and G28 flybys 
occurred when Ganymede was well below the magnetospheric current sheet and the 
magnetospheric magnetic field pointed radially inward towards Jupiter. G29, the last 
flyby of the mission, occurred at a magnetospheric location similar to that previously 
examined on passes G1 and G2. The complete set of passes enables us to examine 
Ganymede in the presence of varying external field orientations. However, some of the 
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passes are not optimal for analysis of internal sources. We will return later in the paper to 
a discussion of criteria used to determine data useful for our purposes. 

The measured magnetic fields contain both the background magnetic field (B ig (t)) of 

Jupiter, which varies very little over the scale of Ganymede’s radius, and the magnetic 
signatures associated with Ganymede, which vary markedly on the same spatial scale. In 
evaluating sources of magnetic field internal to Ganymede, we look for changes 
supplementary to the changes in the Jovian field that would have been measured along 
Galileo’s trajectory had Ganymede been absent. Khurana’s (1997) model of the 
magnetospheric magnetic field provides a rough estimate of the background field, but we 
improve our estimates post-flyby by fitting a polynomial to the components of the 
magnetic field measured before and after the entry into the near vicinity of Ganymede. 

We have assessed the quality of our fit to the background field by using data from 
orbits that do not come close to Ganymede but cut through its orbit at System III 
longitudes similar to those used for this study of the internal field. We fit a polynomial to 
data taken before and after an interval of one hour near the crossing of the Ganymede 
orbit. We find that the polynomial represents the magnetic field measured within the 
excised interval with a standard deviation of less than 2 nT. Furthermore the background 
field changes little over the intervals fitted (near Ganymede’s orbit, the components 
change by <3 nT over distances ~1 Rq= 2634 km) and the variation is roughly linear. 
Errors in the representation of the background field will very slightly affect external 
moments in the fits but internal moments will be negligibly affected. The difference 
between the data and the fitted background field is dominated by the magnetic field from 
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sources within Ganymede and by the effects of magnetopause currents and other currents 
arising from the interaction of magnetospheric plasma with Ganymede. 

Figures showing the measured magnetic field components and the field magnitude 
near closest approach for the first four passes by Ganymede can be found elsewhere 
(Kivelson et al. 1998). Here we plot data from the G28 flyby in order to illustrate the type 
of magnetic field measurements that are used in the analysis. Figure 3 shows the three 
components of the magnetic field and its magnitude from measurements at 0.3 s 
resolution and the fit to the slowly varying background field for 30 minutes around 
closest approach for the G28 pass of May 20, 2000. A schematic view of the field is 
shown in Figure 4 to help put the measurements in context. In Figure 4, the field lines 
are calculated for a vacuum superposition model of the background field of Jupiter 


Fig. 3 


Fig. 4 


(approximated as locally uniform in an orientation appropriate for the G28 pass with 
components given in Table I) and the internal magnetic moment of Ganymede projected 
into the Ganymede-centered plane perpendicular to the flow direction. The G28 
trajectory has been projected into the same plane and marked with time tags. As in 
previously published schematics of the magnetic field (Kivelson et al., 1998), three types 
of field lines can be identified. Those that close on Ganymede at both ends are referred 
to as closed field lines. Those with one end on Ganymede are referred to as open field 
lines, and those with both ends at Jupiter are referred to as Jovian field lines. A 
separatrix between field lines with at least one end on Ganymede and the Jovian field 
lines can be clearly identified. Although in our simplified model no current flows on this 
separatrix boundary, it resembles in other respects the magnetopause. The times of 
magnetopause crossings should be close to the times of the crossings of the separatrix 
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(Kivelson et al. 1997; Williams et al. 1997a, 1997b). However, the cross section of the 
separatrix in the x-y plane is roughly circular and the trajectory lies on a chord, not a 
diameter, of this circle. Thus, Galileo is expected to encounter the separatrix later on its 
inbound crossing and earlier on its outbound crossing than the times for separatrix 
crossings in Figure 4. The times for the crossings are: inbound crossings at -1004:30 
UT in Figure 3 compared with -0959 UT in Figure 4 and outbound crossings at -1019 
UT in Figure 3 compared with —1026 UT in Figure 4, consistent with expectations. 

The separatrix corresponds to the approximate position of the magnetopause but gives 
no insight into the currents that flow on it. However, field rotations evident in the data 
arise because of boundary currents, so their influence can be readily estimated from the 
actual data. Magnetopause currents are sheet-like and produce perturbations that vary 
slowly with distance normal to the boundary. This effect can be seen in Figure 3 where 
an abrupt field rotation at 1004 UT decreases the B y component by 124 nT with a smaller 
change in the B z component and what appears to be a return to background in the B x 
component. (In a minimum variance coordinate system, the change in the maximum 
variance direction across the boundary is 174 nT.) For the next 2 minutes, the B y and B x 
components change very slowly because perturbations arising from magnetopause 
currents dominate those arising from sources internal to Ganymede. Subsequently, the 
field changes rather rapidly as Galileo moves closer to the moon. Analogous signatures 
are familiar for inbound passes through the dayside of the terrestrial magnetosphere 
where the field jumps to roughly double the dipole field magnitude near the 
magnetopause and changes little over several Re until it approaches the dipole field value 
near 6 Re. Inward of this distance, the field magnitude follows the dipole field model 
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rather well. (Figure 1.12 in Russell (1995) shows that a similar pattern is observed on 
passes through the terrestrial dayside magnetosphere.) 

In order to focus on the portion of the pass that yields information useful for 
characterizing internal field sources, we fit only that part of the flyby data in which the 
deviation from background of the field component with the largest change across the 
magnetopause is at least double its jump in magnetic field strength at the magnetopause. 
This approach eliminates portions of the data in which fluctuations arising from the 
proximity of magnetopause currents may contribute extraneous signals while retaining a 
major portion of the signal dominated by sources internal to Ganymede. The relevant 
time intervals are given in Table II. Extending or reducing the intervals by 30 s at each 
end produces less than -10% changes of fit parameters. We fit the data separately for 
each pass with a model containing a centered dipole moment at Ganymede and a uniform 
magnetic field (UFX, UFY, UFZ). We report the magnitude of the magnetic dipole 
moment (in nT) in terms of the magnitude of the magnetic field it produces at 
Ganymede’s surface (i.e., at the radial distance 1 Re) at the dipole equator. (If M is the 
equatorial surface field magnitude, the magnitude of the dipole moment in standard SI 
units of Amp-m 2 is 47 tM(nT)RG(km ) 3 //4 = 1.83xl0 17 A/(nT).) The uniform field 
components approximate the perturbations produced by magnetopause currents, which 
are taken as constant for each separate pass and assumed to vary little over the selected 


Table II 


intervals. The parameters of the internal dipole obtained from the different flybys vary 
considerably with the fits to passes G8 and G7 as the outliers, a matter that we explain in 
the next section. 
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Ganymede’s internal sources fitted with fixed dipole and quadrupole moments 

The dipole components in Table II are inconsistent from one pass to another. As the 
passes occur at different planetocentric latitudes and longitudes, the scattered values may 
reflect contributions of higher order internal multipole moments. Fits of the limited data 
available to dipole plus quadrupole moments requires that the information content of the 
passes be adequate to distinguish the variations related to different multipole moments 
along the trajectory. In order to avoid fitting data to passes with insufficient information 
content, we next assess the sensitivity of the measurements along each of the trajectories 
to contributions of low order multipole moments by assuming that the amplitudes of all 
dipole and quadrupole coefficients are identical at the surface. Assigning a nominal 
surface amplitude of 50 nT to each of the first 8 multipoles ( g,°, gj ,h\,g° 2 ,g\,g\,h\, h \ ) 
as defined in Walker and Russell (1995), we plot (Figure 5) the perturbations that would 
be present along the five trajectories. Most striking is the small amplitude of distinctive 
variations along the G7, G8 and G29 trajectories, all of which have closest approach 
altitudes of 1600 km or more, compared with the signatures on the other, lower altitude, 
trajectories. For the low altitude passes (Gl, G2, and G28), the signatures of some of the 
components are clearly evident, with amplitudes of at least 20 nT. Therefore we identify 
the latter three passes as relevant to the determination of the internal dipole plus 
quadrupole coefficients. In addition, each pass includes contributions from magnetopause 
currents that we approximate as uniform fields (which appear in the multipole fit as the 
first order external coefficients: G,°, G\,H \ ). 

The values of the uniform field contributions are expected to vary from pass to pass, 
but the multipole moments are assumed not to vary over time scales of months and years. 


Fig. 5 
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A modification of the standard least squares fitting technique described in Appendix A 
allows us to combine data from all three flybys as a single data set to determine the 
internal moments. In an unweighted least square minimization, this approach will 
optimize the fit to the data from passes that recorded the largest field magnitudes. Some 
accommodation is required to account for the fact that the flybys occurred at different 
altitudes and that the maximum field strength encountered changed markedly from one 
pass to the next. In order to balance contributions from different passes, we weight data 

from pass i inversely with B /(^.(5^ x ) 2 j , where B‘ mK is the maximum field strength 


in the ith interval (Table II). 

Broad coverage in latitude and longitude are needed to distinguish different multipole 
components. G1 and G2 came closest to Ganymede at similar planetocentric longitudes, 
whereas closest approach on G28 was separated by about 150° in longitude (Table I and 
Figure 1). Thus G28 data play a unique role in the specification of the longitudinal 
structure of the field. By double-weighting the G28 data, thus partially compensating for 
the imbalance produced by having two passes at roughly the same planetocentric 
longitude, we improve our ability to distinguish among various multipole coefficients. 

In the multipole fit, we model the data set with 17 parameters: 3 different sets of 3 
parameters to represent the uniform fields that apply to the individual flyby portions of 
the data set and 8 additional parameters that characterize the dipole and quadrupole 
moments. Values of the multipole moments that we obtain are given in Table III. 

Various features of this fit are worthy of comment. The dipole moment (magnitude 
713 nT, 176° colatitude, with southern hemisphere pole rotated 25° from the Jupiter- 
facing meridian) remains close to that found in early estimates. The uniform field 


Table III 
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contributions are readily understood in relation to aspects of plasma currents in the 
interaction region. In particular, the Alfven wing bendback accounts for negative B x 
perturbations above Ganymede’s equator (passes G1 and G2) and positive B x 
perturbations below the equator (pass G28). The B z component should be little changed 
across the magnetopause in high latitude passes such as G1 and G2, as evident from 
Figure 4 and consistent with small UFZ as in Table III, but requires field rotation across 
the magnetopause for low latitude passes such as G28 where Table III indicates that UFZ 
= 80 nT. 

Particularly notable is the very small ratio (0.0016) of the power in the second order 
internal multipole moments to that in the first order moments. The ratio of power in the 
quadrupole to the power in the dipole coefficients has been used elsewhere (Elphic and 
Russell 1978; Connemey 1993) to infer the depth of the sources of magnetic field 
assuming dynamo field generation. This is a subject to which we will return in the 
discussion. 

Ganymede’s internal sources fitted with a fixed plus induced dipole moment 

Evidence of inductive magnetic fields has been found at Europa and Callisto (Khurana 
et al. 1998; Kivelson et al. 1999, 2000; Zimmer et al. 2000). The cited references 
explain the observed magnetic signature in terms of the inductive response of an 
electrically conducting shell to a time varying magnetic field. Within Jupiter’s 
magnetosphere, the orientation of the external magnetic field present at a moon varies at 
Jupiter’s synodic period as viewed from that moon. The component along the direction 
radial towards Jupiter exhibits the largest variation. If an electrically conducting shell is 
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present within the moon, the time varying field drives inductive currents within that shell. 
These induced currents produce a time varying magnetic moment that lies in the moon’s 
spin equatorial plane and points roughly towards and away from Jupiter. 

The induced magnetic field signatures at Europa and Callisto are very clear. A similar 
signature near Ganymede will be hidden in the signature created by the internal dipole as 
noted. Figure 6 and Table I show that the G28 pass occurred when Ganymede was 
located well below the center of Jupiter’s magnetospheric plasma sheet where the 
externally imposed magnetic field tilted radially inward towards Jupiter. G1 and G2 
occurred at times when the externally imposed magnetic field tilted radially outward from 
Jupiter. Thus, the three flybys used to obtain low order multipole coefficients are also 
valuable in determining whether an inductive response is present. The external field for 
G1 and G2 was directed radially outward, whereas the external field for G28 was directed 
radially inward. Were there an induced magnetic moment, it would be antiparallel to the 
radial (relative to Jupiter) component of the exiemai field and its orientation for G28 
would differ from that for G1 and G2. 

We use a two step approach to test the possibility that the variations among the fits to 
the three critical passes given in Table II arise through inductive responses. First we 
assume that the dipole moment does not vary. We use the approach outlined in the 
appendix but this time to find a single best-fit dipole moment plus different uniform 
fields for each orbit. Weighting by the inverse of the maximum field measured in each 
interval and doubling the weight of the G28 interval as before, we follow the procedure 
previously summarized but this time determine only 12 fit parameters. Results for this fit 
are given in Table IV which shows changes of 2% in Mz , 17% in Mx and 13% in My 


Fig. 6 
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relative to the values of Table III. The total rms error (see Appendix C for definition) 
increases from 13.5 nT to 15.1 nT. As anticipated, the 3 dipole parameters fit the data 
less accurately than do 3 dipole and 5 quadrupole parameters. I 

We next modify our approach by assuming that the internal sources include both a 
time varying moment arising from induction and a constant dipole moment. As described 
above, a time-varying uniform magnetic field imposed on a perfectly conducting sphere 
induces a magnetic dipole with equatorial surface magnitude half that of the driving field, 
oriented antiparallel to the instantaneous driving field. Imperfect electrical conductivity 
decreases the magnitude of the induced dipole moment and introduces a phase lag. 
Zimmer et al. (2000) have applied the theory of conducting spherical shells to the 
problem of the inductive responses of Europa and Callisto. They show that a shell of 
Callisto’s radius (which is close to Ganymede’s) responds with a phase lag corresponding 
to less than 10° of Jupiter rotation if the shell thickness is more than ~10 km, provided its 
conductivity is greater than or equal to that of terrestrial seawater. Such a phase shift is 
too small to be detected, particularly for encounters that occur well off the center of the 
magnetospheric current sheet in a region where the field changes slowly with time. 
Consequently, we carry out our analysis ignoring corrections for phase lag and 
characterize the inductive response in terms of the amplitude of the induced field. We 
introduce a parameter, the response efficiency a , that equals the ratio of the magnitude 
of the induced dipole moment at Ganymede’s surface to that of a perfectly conducting 
sphere of 1 Rq. As the dominant time variability of the background field is in the 
direction radial from Jupiter (the h-di recti on in GphiO coordinates), we assume that Mx 
and M z are fixed for all passes and that My(t) responds to the driving field, B b Y s {t) as 
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My(t) = M Yo - aBy g (t ) . Here M Vo is the T-component of the constant dipole moment 
and the inductive contribution is opposed to the driving field, which we take to have the 
values provided in Table I. Values of a < 1 can arise from a combination of imperfect 
electrical conductivity, a small shell thickness, and/or a shell radius <1 Rg- 

With the introduction of the response efficiency, at least four fixed internal parameters 
must be specified, the three components of the permanent magnetic moment (where 
M Yo is g' t of the dipole fit) and a. Again we use the approach described in Appendix A, 
weighting all passes with the inverse of the maximum field in the interval and double 
weighting the G28 pass. Table V gives the results of this analysis. 

The results show that M Yo obtained in this final analysis, is slightly smaller than the 


value obtained by using the three passes to infer the components of the dipole field 
(Table IV) and close to the value of My found by fitting both first and second order 

multipoles (Table III). 

The rms error of fit to a dipole with varying My is 11.5 nT, whereas the rms error is 
13.5 nT in the quadruple fit (Table III, footnote 4). (In carrying out the least squares fit, 
we use weighted rms errors, but in the summary tables we provide the rms errors without 


weighting and in Table VI we give both.) The improvement in the rms values is 
significant because the model based on an inductive response requires only 4 internal 
parameters whereas the quadrupole fit uses 8. 

In Figure 7 we plot M r (t) = M Yo - aB b /(t ) for the inferred response of a = 0.84 as a 

function of the induced moment that would be found for an inductive response with a - 1 
which we refer to as Modeled My. Here My 0 is the permanent moment. The dark line 
shows the expected values for a = 0.84, the value used to place the data points, and the 


Table VI 


Fig. 7 
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gray line gives the values for a = 1. Although such a large value of a lies well outside 
the uncertainty allowed by the formal error analysis (±0.02), we believe that the formal 
error analysis underestimates the uncertainty of the fit parameter. In the next section we 
discuss why we believe that a more realistic assessment of the uncertainty arising from 
the limited number of passes and the inadequacies of the model would yield larger 
uncertainties of the fit parameters. We cannot, therefore, exclude the possibility of a as 
large as 1 . 

Discussion 

We start by considering whether our analysis provides conclusive evidence that 
Ganymede’s field is the sum of a permanent magnetic moment plus and an inductive 
moment. We remain cautious about such a conclusion as the data are also well fitted by 
the eight-parameter dipole plus quadrupole model. The quality of the fits is discussed in 
detail in Appendix C where relevant quantities such as condition number and weighted 
and unweighted rms error are defined. If summed over the fitted passes, the model of 
Table V gives an rms error of 1 1.5 nT and a weighted rms error of 13.2 while the model 
of Table III give an rms error of 13.5 nT and a weighted rms error of 12.6. These values 
are sufficiently close to one another that we cannot conclude that one model is to be 
favored over the other. Nonetheless, we believe that the success of a four-parameter fit 
based on an a priori hypothesis must be more than coincidence. The rms errors of dipole 
moments for individual passes in Table V (fixed plus induced dipole moment with an 
84% inductive response efficiency) is only slightly greater than the pass by pass rms 
errors in Table II (dipole moments determined independently for each pass), and the 
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considerable variation among the moments obtained in Table II yields the improbable 
result that a nominally permanent dipole moment changes between passes or that the 
errors are comparable with the differences between the estimates. The latter interpretation 
is inconsistent with errors as large as ten times the formal error estimates. 

Although we used the digitization step size of 8 nT as the basis for estimating errors in 
the tables (see discussion in Appendix C), the actual measurement errors averaged over 
spin periods are significantly smaller. The total rms errors of about 12 nT over the 
different passes used for fitting the permanent moments for both the inductive and 
quadrupolar models exceed the measurement errors and arise from other effects. 
Contributing to the rms errors are higher order internal moments of Ganymede, 
fluctuations produced by local plasma currents and temporal variations of the Jovian 
magnetosphere on scales of minutes. Indeed, we suggest that because of the large number 
of extraneous contributions to the measurements, the formal errors in the fit parameters 
underestimate their actual uncertainty. 

Error analysis, then, is consistent with either the dipole plus induced field model or the 
dipole plus quadruple model. In the quadruple fit, only two of the quadrupole 
coefficients are larger than the uncertainty in their values. If a theoretical argument were 
to be developed that singled out the quadrupolar coefficients g 2 1 and h 2 2 as the principal 
higher order moments, one could look to Table III for support, but we are not aware of 
any arguments that other moments of the same order should be unimportant. But, we 
reiterate that we favor the inductive response model because we find it hard to suppose 
that the dominant quadrupolar moments fortuitously produce perturbations on the three 
orbits that match those of a time varying induced response. 
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If one accepts the proposal that there is an inductive response, one needs to consider 
where the induced currents are flowing and why the response may be at a level lower 
than "ts theoretical upper limit. As for Europa (Kivelson et al. 2000; Zimmer et al. 2000), 
the icy outer layer of Ganymede cannot provide sufficient electrical conductivity to 
account for the response unless there is a melted layer, a water ocean, below the surface. 

If a water layer is present below but near the surface, an 84% response is consistent 
with a range of models of the conducting layer. Assuming the conductivity of terrestrial 
sea water, the upper limit amplitude is consistent with a conducting layer of greater than 
~1 km thickness very near the surface (Zimmer et al. 2000), but because the responses 
would be reduced in amplitude for some of the passes if phase lag were considered, the 
layer would have to be closer to 1 0 km in thickness. Alternatively, a thick conductive 
layer may be buried more deeply beneath the surface, the interpretation that we favor. 
Taking into account the expected cubic falloff with distance from the conducting surface, 
a thick layer buried at depth (1-(0.84) 1/3 ) = ~150 km would account for the reduction 

of the magnetic signal to 84% of its maximum possible value. 

Schubert et al. (1986), invoking the equivalence between pressure and depth below 
Ganymede’s surface, converted the melting temperature vs. pressure curve (the liquidus) 
of pure water ice to a curve displaying melting temperature vs. depth beneath 
Ganymede’s surface. Their plot forms the basis of Figure 8, which shows that in the 
interior of Ganymede the minimum temperature at which water ice melts occurs at a 
depth of 170 km. (The minimum melting temperature would decrease if salt were 
dissolved in the water.) We have added a possible temperature profile to the plot of the 
liquidus. If the temperature rises rapidly in the first 100-150 km beneath the surface and 


Fig. 8 
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increases more slowly (adiabatically) at greater depth, a layer of melt would be found 
sandwiched between layers of ice surrounding the -170 km depth. Thus the most 
probable value of the response factor a suggests a depth for the buried ocean that is 
reasonable on physical grounds. Additional work is underway to relate the observations 
to properties of a possible subsurface ocean. 

Finally we return to the implications of the small values obtained for the quadrupole 
coefficients (Table III). (Note that the quadrupole moments are smaller if one also allows 
for inductive responses, but the values obtained for a fixed dipole moment will suffice to 
make the arguments.) One possible interpretation of the small amplitude of the higher 
order multipoles is that the region in which the permanent internal field is generated is 
deeply buried. This concept, proposed by Lowes (1974), has been applied by Elphic and 
Russell (1978) to a determination of the “apparent source depths” of the fields of 
Mercury and Jupiter and by Connemey (1993) to Earth and the outer planets. The 
approach identifies a depth at which the power in multipoles of order n 
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is comparable for all orders, based on the idea that the moments arise from turbulent 
motions in the source region. Although the ratio of first two orders of multipole 
moments can be no more than suggestive, we note that the ratio of power is 
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which yields 0.0016 from Table III, down by almost two orders of magnitude relative to 
the analogous ratios for Earth or Jupiter (Elphic and Russell 1978; Connemey 1993). 
With r as defined as the radius of the source sphere, i.e., the distance from Ganymede’s 
center to the “apparent source depth” at which U 2 / U\ = 1, 


La 

R, 


_m~0 

2 £[(*,' 

V rn =0 J 


v 1 / 2 


( 3 ) 


For the values in Table III, r ns = 0.049 Rq =130 km. This estimate of the apparent source 
radius is hard to reconcile with estimates of the size of the metallic core whose radius is 
believed to fall between 0.25 Rg — 658 km and 0.5 Rq = 1317 km (Anderson et al. 1 996). 
It should be noted, however, that Elphic and Russell (1978) find that the quadrupole 
power at Earth is lower by almost an order of magnitude than the trend provided by the 
dipole moment and the n > 2 multipole moments, so the estimate of the source depth 
from U\ and Uj alone may be misleading. It does seem reasonable to conclude, however, 
that if Ganymede’s field is generated by dynamo action, the source lies far below the 
surface. 


The small ratio of the quadrupole to dipole power does not necessarily require a 
deeply buried dynamo because dynamo models admit a great variety of solutions. 
Magnetoconvection, for example, provides an interesting alternative to the dynamo 
mechanism for generating the fields of satellites (e.g., Sarson et al. 1997). 
Magnetoconvection produces a planetary field by amplification of an externally imposed 
magnetic field such as the average at Ganymede. Fluid motions amplify the imposed 
field much as in a dynamo, but the internal field decays if the external field vanishes. A 
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possible consequence of this model could be that the symmetry of a uniform external 
field imposes a dominantly dipolar symmetry on the internal field and therefore produces 
an internal response with weak higher order multipoles. Yet another proposed source of 
the permanent magnetic moment is remanent ferromagnetism (Crary & Bagenal 1998) 
which would also imply dipolar symmetry for the field. Thus our conclusion that the 
quadrupolar terms are very small may be regarded as providing encouragement for 
further investigation of magnetoconvection or remanent magnetization as the mechanism 
for generating Ganymede’s permanent magnetic field. 

Conclusions 

Magnetometer data from all five passes by Ganymede have been used to refine our 
evaluation of the moon’s permanent dipole moment. The fits of Tables III and V are 
consistent with a magnetic moment characterized by an equatorial surface field 
magnitude M = 719 ± 2 nT tilted by 176°± 1° from the spin axis with the pole in the 
southern hemisphere rotated by 24°±1° from the Jupiter-facing meridian plane towards 
the trailing hemisphere. In SI units, the dipole moment is = 1.32xl0 20 Amp-m 2 . The 
magnitude is very slightly smaller than initially suggested (Kivelson et al. 1996) and the 
orientation is closer to antiparallel to the spin axis. 

Magnetometer data from Galileo’s multiple flybys of Ganymede provide significant 
but not unambiguous, evidence that the moon, like its neighboring satellites Europa and 
Callisto, responds inductively to Jupiter’s time varying magnetic field. The response is 
roughly 84% of the response expected for a metallic sphere of the same radius but there 
may be considerable uncertainty in this estimate. The outer layer of Ganymede is 
predominantly formed of water ice, which is not a good conductor of electricity. The 
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most likely source of electrical conductivity within a water ice environment is a liquid 
water layer bearing electrolytes such as salts and acids. Our data indicate that the locus 
of the current carrying layer is buried within the 800 !km thick outer icy shell (Anderson 
et ah 1996). Melting within the ice shell occurs most readily near 170 km beneath the 
surface where increasing pressure reduces the melting temperature to its lowest value. If 
we assume conductivity comparable to or higher than that of terrestrial sea water, our 
analysis is consistent with a depth for the conducting layer close to this location. The 
ocean would lie between two layers of ice. Ganymede’s large size relative to its two ice 
covered neighbors Europa and Callisto and the evidence previously found that the nearby 
moons are responding inductively lends credibility to the evidence that Ganymede is 
responding inductively to the temporal variations of the Jovian magnetic field. 

Although we cannot rule out other interpretations of the data discussed above, we 
view the inductive response interpretation as highly probable. Regrettably, as illustrated 
in Figure 5, data from G29, the final Ganymede encounter of the Galileo mission, did not 
add meaningfully to the information already available for the analysis of Ganymede’s 
internal magnetic properties. However, planning for more Jupiter orbiting spacecraft 
missions continues, and it should be recognized that future opportunities to take 
magnetometer data in Ganymede’s environment can be designed to acquire data 
necessary to distinguish between the signature of fixed internal quadrupole moments and 
the signature of an induced dipole moment. 
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Appendix A: Least squares fits 

We use a weighted least squares fit to determine a single set of multipole components 
for multiple passes and d ifferent constant field components for each individual flyby. The 
input data are values of the magnetic field b along each of a set of trajectories. These 
measurements are related to the matrix of the model field coefficients x by a matrix A, a 
function of Galileo’s location relative to Ganymede that varies from one pass to another. 
With the error defined as: e = (Ax - b), a least squares fit requires us to minimize the 
square sum of E = e r e/2 (where the superscript T indicates the transpose) with respect to 
the unknown parameters of the model. This requires dE/dx = 0. The parameters of the 
model are then the solutions of 

^(Ax-b) r (Ax-b)=^A r (Ax-b) + ^(Ax-b) r A = A r (Ax-b) = 0, (Al) 
which by matrix inversion gives: 

x = (A r A)~‘ A r b . (A2) 

As the passes vary considerably in altitude and in maximum field strength, fits to the 
full data set are biased towards fitting the lowest altitude passes (see Table I). In order to 
optimize the fit over the full data set, we weight the contributions of different passes 

inversely by B‘ mx ) 2 )' , where is the maximum field strength in the /th 

interval. In our application, the contribution of G28 is treated as if there had been two 
identical G28 passes, so i = 1 ,.. .4. In a weighted least squares fit one redefines E as: E = 
e W e/2, where W is a diagonal weighting matrix. The equivalent to (A2) is then: 

x = (WA T WA)-' WA T Wb . ( a3) 
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In this paper the measurements b consist of the three components of the magnetic field 
in the GphiO coordinate system. If we fit only dipole coefficients of the internal field, the 
matrix A for each data point is a 3x12 matrix, built up in the following form: 
[ 3 A, UF(G7), UF(G2), UF(G7), UF(G2S)] where 3 A is a 3x3 submatrix that relates the 
dipole field along the orbit to the Cartesian components of the dipole moment: 

R = yfx T +Y T +Z T 

3 A u = (3X 2 -R 2 )/ R 5 , 3 A n = 3XY / R 5 , 3 A x3 = 3XZ/ R 5 , (M) 

3 A 2] = 3XY / R 5 , 3 A 22 =QY 2 -R 2 )/R 5 , 3 A n = 3YZ/ R 5 , 

3 A 31 = 3XZ/ R 5 , 3 A 32 - 3YZ / R 5 , 3 A 33 =(3 Z 2 -R 2 )/R 5 , 

The next three 3x3 submatrices are I if the data being fitted correspond to the 
specified orbit, or 0 if not. 

For each data point this matrix is calculated and added as new rows, which, for n data 
points leads to a new matrix A of dimension 3nxl6. The above calculation of the least 
squares fit is then performed as specified by equation (A2). The extension to fits of 
higher order multipoles produces larger matrices, but there is no conceptual change in the 

approach. 

For the final fit, a parameter of a different type (a) is fitted in addition to the dipole 
coefficients but the process remains as described above. 

Appendix B. Singular Value Decomposition 

We have used the generalized inverse theory to reconfirm the least squares analysis and 
to determine the uniqueness and the robustness of the fits. The theory of generalized 
inversion was discussed in detail by Lanczos (1961) and highlighted by Jackson (1972) in 
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the context of geophysical data inversion. Applications of the technique to invert 

geophysical magnetic data have been made by Pedersen (1975) and Connemey (1981). 

I 

The linear set of equations that we must invert can be written in the matrix form as: 

y = ax (Bl) 

where Y is a column matrix of the N magnetic field observations, A is the NxM matrix 
which relates observations to the model parameters X (the various internal and external 
spherical harmonic terms and the response). (In terms of the index n used in Appendix 
A,N=3n.) As N >M in our case, the linear system is overdetermined and a least squares 
inversion provides the best estimates of the model parameters. 

The idea behind the generalized inverse technique is to decompose the A matrix into 
the following form: 


A = U A V T 

(B2) 

where U, A, and V satisfy the following eigenvalue problems: 


aa t u = A 2 U 

(B3) 

A t AV = A 2 V 

(B4) 


The A matrix is a diagonal MxM matrix with the M eigenvalues arranged along the 
diagonal of the matrix in order of decreasing magnitude. U is an NxM matrix with the M 
eigenvectors of equation (3) arranged in //columns; V is an MxM matrix with the M 

eigenvectors of equation (4) arranged in M columns. The least squares inverse of matrix 
A is 

H = VAU T (B5) 
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Appendix C. Estimating the Validity of the Fits: Condition Number and RMS 
Deviations 

An examination of the U, A, and V matrices gives insight into the inversion process. 
For example, the ratio of the largest to the lowest eigenvalue (called the condition 
number of the matrix A) provides a measure of the invertibility of matrix A. If the 
condition number is large, the solution is poorly constrained by the data and the poorly 
determined model parameters have large errors associated with them. There is no 
rigorous selection criterion for how large the condition number can become and still 
imply a good fit. Indeed the number is expected to increase with the number of 
parameters being determined, yet this does not mean that the solutions are poor. In 
applications to internal fields of planets, values as large as 60 are regarded as adequately 
small. The solutions of Tables III-V have condition numbers smaller than 32.6 and are 
tabulated in Table VI. It can be seen that for least square fits in which quadruple terms 
are omitted (the first two rows of the table), the condition numbers are quite small (< 13). 
The condition number increases to a value of 32.6 when dipole, quadrupole and external 
field parameters are optimized. When the induction response is simultaneously 
optimized, the condition number almost doubles, increasing to a value of 51.4. This large 
increase reflects the fact that the available data cannot distinguish the h\ quadrupole 
contributions from the induction response. The V matrix for this fit reveals that these two 
parameters require large contributions from the three eigenvectors corresponding to the 
three smallest eigenvalues. 
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Estimates of the errors of the fit parameters can be obtained from the matrix V. The 
standard error associated with a parameter j (assuming that the observations are 
statistically independent) is given by: 


Sj — <j B 




(Cl) 


where cr B is the standard uncertainty associated with the measurements (assumed 
constant in this work). The uncertainties of the fit parameters given in Table III- V were 
calculated from equation (B6) for P = M and a B = 8 nT. This value of the standard 
uncertainty corresponds to the digitization of the magnetometer in its highest field range 
[Kivelson et al. 1992] and overestimates the error of measurement. Uncertainties in 
position relative to Ganymede are less than 1 km, implying that the uncertainty in range 
A R/R = 4x10^ is allowed for within the value of 8nT that we used for 

It is useful to confirm that the model provides a good representation of the data by 
calculating the root mean square deviation of the fit to the measured quantities. In Table 
VI we tabulate the rms deviations both for weighted and unweighted cases, where 
weighting is relevant to the approach used to reduce the disproportionately large 
contribution of the lowest altitude flyby. The definition of the summed weighted rms of 
fit ( rmsw ) is: 

I" 

rmsw - J — J 

X(w, +w 2 +2u' 2g ) l/2 


(C2) 
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where the weight of the ith pass by Ganymede (Gi) is w ( . [5 max (GOr and B max (GO is 

tabulated in Table II. Here A B jGi represents the deviations of the Not individual data 
points on pass Gi from the fit for that pass. N = N G ,+ N G2 + 2 N G28 is the total number 
of data points in the data set. The summed rms of fit is defined by 

rms = rmsw with vv / = 1 for all i (C3) 

Lanczos (1961) has shown that a more robust inverse (called the generalized inverse) 
is obtained by retaining only the P largest eigenvalues of the V matrix and their 
eigenvectors in equation (B5). However, in some cases, the eigenvector corresponding to 
the smallest eigenvalues contributes significantly to the fit and it may not be appropriate 
to drop it. For the Ganymede fits, we found that when we dropped the smallest 
eigenvalue from the 18 parameter solution (fitted parameters listed in the last row of table 
VI), the new solution was almost identical to the 17 parameter solution of the third row. 
The summed rms error of fit changed from 13.5 to 13.4. The insignificant change 
confirms the previously noted fact that the inductive response merely reproduces 
contributions that could equally be attributed to h\. In other cases, when the smallest 
eigenvalues was dropped, the rms error increased by almost a factor of 2. For example, 
the rms error of fit to the case of Table III increased from 13.5 nT to 28.7 nT when the 
eigenvector with the smallest eigenvalue was dropped, indicating that it was not 
appropriate to use the generalized inverse solution. 
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Table I; Information regarding Galileo’s encounters with Ganymede. 

pass date C/AUT LT Alt.(km) | Planetocentric U) Location 2 Background field 3 

lat.n long. (°) rel.tops B^x 

Gl 06/27/96 0629:07 11.27 838 30 247 above 6 -79 -79 

"G2 09/06/96 1859:34 10.78 264 79 236 above 17 -73 -85 
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The rms error is calculated for the difference between data and model in the fitting interval. 

Arguments in the text explain that the information in the G7 and G8 passes does not constrain the dipole moment, so these passes are omitted from the averages in Tables III-V. 









Table III: Fits to dipole and quadrupole moments fixed for all passes, 
with uniform fields varying from pass to pass . 1 
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Table V: Fits to a dipole moment and induction efficiency factor a fixed for all passes 
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Figure 1. Plots of Galileo’s passes by Ganymede. The coordinate system is defined in 
the text. 

Figure 2. Schematic of the coordinate systems used in the analysis. The origin of the 
axes is at the center of Ganymede. 

Figure 3. Magnetometer data for G28 (three components and the field magnitude in nT) 
vs. UT. The rotations on a time scale of a few minutes near 1005 and 1019 UT (see 
shading) represent crossings of the current layer of Ganymede’s magnetopause. 

Figure 4. Vacuum superposition of an internal dipole field and a uniform field with 
orientation appropriate to that present during the G28 pass. The cut is in the yz plane of 
the GphiO coordinate system and Galileo’s trajectory has been projected into this plane. 
Following the nomenclature explained in the text, Jovian field lines are dashed, closed 
field lines are solid, open field lines are dotted, and the separatrix is a heavy solid curve. 
Figure 5. Traces of the contributions of the dipole and quadrupolar moments along the 
portions of the passes used for fitting internal moments assuming that the surface 
amplitudes of all coefficients at the surface are 50 nT. G1 is blue, G2 is red, G7 is green, 
G8 is cyan, G28 is black, and G29 is gray. 

Figure 6. The variation of modeled field components transverse to Ganymede’s spin axis 
over a Jupiter rotation period in Ganymede’s frame (Khurana 1997). The values 
anticipated at the times of the Ganymede passes are indicated. 

Figure 7. Inferred My vs. the upper limit My that can be induced in a spherical conductor 
of radius IRq labeled as Modeled My. A dark straight line represents the response to 
induction assuming a permanent My = 49 nT and a response of 82% of the response of a 
metallic sphere of radius 1 Rg to By imposed on Ganymede by the time- varying Jovian 
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magnetospheric magnetic field. The gray line represents the response for a metallic 
sphere. Points corresponding to the passes used in the analysis are shown. 

Figure 8. Plot of the melting curve (liquidus) of pure water ice and the boundaries 
between different forms of ice as a function of depth beneath the surface of Ganymede 
(after Schubert et al. 1986). A temperature profile as a function of depth beneath the 
surface consistent with having a region of melt near 170 km depth has been 
superimposed. 
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Galileo Ganymede Flyby Trajectories (GphiO) 
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to Jupiter 


Figure 2 
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